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ABSTRACT 

We present HST/WFPC2 observations across the disk of the nearby isolated dwarf SO galaxy 
NGC 404, which hosts an extended gas disk. The locations of our fields contain a roughly equal mixture 
of bulge and disk stars. All of our resolved stellar photometry reaches tufsuw = 26 {Mfsiaw — —1-4), 
which covers 2.5 magnitudes of the red giant branch and main sequence stars with ages <300 Myr. 
Our deepest field reaches tofsi41V = 27.2 (Mfsuw = —0.2), sufficient to resolve the red clump and 
main-sequence stars with ages <500 Myr. Although we detect trace amounts of star formation at 
times more recent than 10 Gyr for all fields, the proportion of red giant stars to asymptotic giants and 
main sequence stars suggests that the disk is dominated by an ancient (>10 Gyr) population. Detailed 
modeling of the color-magnitude diagram suggests that ~70% of the stellar mass in the NGC 404 disk 
formed by z ~ 2 (10 Gyr ago) and at least ^90% formed prior to z ~ 1 (8 Gyr ago). These results 
indicate that the stellar populations of the NGC 404 disk are on average significantly older than those 
of other nearby disk galaxies, suggesting that early and late type disks may have different long-term 
evolutionary histories, not simply differences in their recent star formation rates. Comparisons of the 
spatial distribution of the young stellar mass and FUV emission in GALEX images show that the 
brightest FUV regions contain the youngest stars, but that some young stars (<160 Myr) lie outside 
of these regions. FUV luminosity appears to be strongly affected by both age and stellar mass within 
individual regions. Finally, we use our measurements to infer the relationship between the star for- 
mation rate and the gas density of the disk at previous epochs. We find that most of the history of 
the NGC 404 disk is consistent with star formation that has decreased with the gas density according 
to the Schmidt law. However, ~0.5-l Gyr ago, the star formation rate was unusually low for the 
inferred gas density, consistent with the possibility that there was a gas accretion event that reignited 
star formation ~0.5 Gyr ago. Such an event could explain why this SO galaxy hosts an extended gas 
disk. 

Subject headings: galaxies: individual (NGC-404) — galaxies: stellar populations — galaxies: spiral 
— galaxies: evolution 



1. INTRODUCTION 

While both redshift surveys and stellar population 
studies agree that late-type disks form the majority 

of th ei r stars by z ~ 1 ( e.e.. iRavindranath et all 

2004 iBarderi et al.l 120051: iMelbourne et all 120071; 
Williams et al. I I2009aft . the ages of early-type disks 
are not as well constrained. The theoretical and high 
redshift observational data suggest that the old popu- 
lations of early type disks may be systematically older 
than those of later type disks. On the other hand, the 
morphological differences between early and late type 
disks may be of recent origin, due only to differences in 
star formation over the past few hundred Myr. Such 
a scenario has has been suggested for the SO galaxy 
NGC 5102 (|Davidgell2008T ). as well as for the differences 
between dEs and dlrrs. 

Unfortunately, little is known about the evolution- 
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ary history of SO galaxies, especially in the field. Al- 
though typical SO galaxies appear to differ significantly 
from E galaxies, it is not yet cl ear where SO galaxies 
fit into the Hubble tuning fork (jvan den Berghl [200l 
or w hat amount of st ar formation they typically con- 
tain dTemi et al.H200cl . While non-interacting SO galax- 
ies tend to have di sks similar in thickness to those o f 
late-type galaxies ()Pohlen et al.l 12004 iLi et al.l [20091) . 
other properties of SO disks differ significantly from 
those of late-type galaxies. For example, their color 
gradie nts are less pronounce d than those of late-type 
disks (jMacArthur et al.ll2004) . an d SO galaxies are typ i- 
cally fainter than later type disks (jvan den Berghl feOOg) . 
SO galaxies are predominantly found in galaxy clus- 
ters and are mostly strip ped of H I ()Poggianti et al.l 
120091 : Ivan den B crgh 2009) , whereas field SO galaxies like 
NGC 404 are les s common and could be a type of post- 
starburst dwarf ()Davidgeil2008[ ). 

As the nearest face-on example of an SO galaxy 
(3.05 Mpc, m-M = 27.42 ± 0.02; M K = -18.9 ; 
A v =0.08; 2=11°; K = 190 km s-^Dalcanton et al.ll2009l; 

I2004J), NGC 404 has 



ISchlegel et all 119981 : Idel Rio et al 



been previously studied in some detail. Although its 
optical morphology suggests that NGC 404 is a classic 
"red and dead" early-type galaxy, the center appears to 
have a significant young co mponent. With long-sli t spec- 
troscopy of the galaxy core. iBouchard et al.l J2010) found 
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evidence for intermediate-aged, and old stellar popula- 
tions, along with a much younger component that fades 
very quickly with radius. The nucleus also h as a domi- 
nant 1 Gyr old population (^50% of the mass, Scth ct al. 
2010). UV spectroscopy also indicate s some very youn g 
(< 10 Myr) stars in the galaxy center (jMaoz et al.lll998l ). 
Farther out, NGC 404 has a very large (>20 kpc) stel- 
lar disk that is know n to be dominated by red giants 
(Tikho nov et al"]|2003t ). This disk also contains a large 
amount of H I gas in a w ell-defined inner dis k (<5 kpc) 
and a warped outer disk (|del Ri'o et al.l[20 04). although 
the density of gas is low (<2x 10 20 M Q pc -2 throughout) 
and only a faint ring of star formation is detected in the 
far ultra-violet (FUV) 5 dWiklind fc Henkellfl990l) . 

The origi n of the gas disk of NGC 404 is somewhat 
mysterious. I del Rio et al.l (|2004T ) suggest that the disk is 
likely due to a recent (0.5-1 Gyr ago) merger event with 
a dlrr galaxy, motivated by their detection of a possible 
warp in the disk. However, NG C 404 is very isolated for a 
galaxy of its mass (0=-l.O; see lKarachentsev et al.ll2004_ 
for details). With no other galaxies of any kind detected 
within 1 Mpc of NGC 404, such a merger is not obvi- 
ously consistent with the galaxy's environment. The gas 
disk could therefore simply be due to late time gas infall 
from filaments rather than accretion of bound objects. 
Furthermore, a gas disk of such low density could have 
a very long lifetime given the low star formation rate, 
suggesting that it could have been accreted far earlier. 
By comparing the stellar populations and gas densities, 
we can constrain the origin of this unusual disk. 

We examine the stellar populations of NGC 404 in 
detail using deep observations from HST. We deter- 
mine the star formation history (SFH) of several por- 
tions of the galaxy by fitting the distribution of stars 
in color-magnitude diagrams (CMDs) with model dis- 
tributions deter mined from st ellar evolution isochrones 
(as described in [Dolphin 2002]). Section 2 describes our 
data set and analysis procedures. Section 3 presents the 
results of our measurements. Section 4 interprets the 
measurements in the context of the gas properties and 
our understanding of SO galaxies, and Section 5 summa- 
rizes our conclusions. W e assume a distance of 3.05 Mpc 
(|Dalcanton et al.l [20091 for conversions of angular mea- 
surements to physical distances and adopt an inclination 
angle i—ll° (|del Rfo et al.|[2004l ) for surfa ce density mea- 
surem ents. We adopt a five-year WMAP (|Dunklev et al.l 
2009) cosmology for all conversions between time and 
redshift. 

2. DATA ACQUISITION, REDUCTION, AND ANALYSIS 

2.1. Acquisition 

From 2007-Aug-08 to 2007-Sep-20, we observed a field 
in the NGC 404 disk located at R.A. (2000) = 17.32325 
(01:09:17.6), decl. (2000) = 35.74856 (+35:44:55) with a 
rotation angle PA_V3=50.0 degrees. From 2009-Feb-16 
to 2009-Feb-20, we performed shallower observations for 
2 fields located at R.A. (2000) =17.36697 (01:09:28.1), 
decl. (2000) = 35.76117 (+35:45:40) with a rotation angle 
PA_V3=230.0 and R.A. (2000) =17.333368 (01:09:20.0), 
decl. (2000) = 35.70205 (+35:42:07) with a rotation an- 
gle PA_V3=230.0. Figure Q] shows outlines of the fields' 

5 http:/ /www. galex.caltech.edu/newsroom/glx2008-02r. html 



locations. Our field locations were chosen to maximize 
the number of disk stars and avoid crowding. 

In the deep fiel d, we obtained 15 full-orbit exposures 
with the WFPC2 (jFord et al.ll 19981) through the F606W 
(wide V) filter, and 29 full-orbit exposures through the 
F814W (/ equivalent) filter. These data totaled 39000 s 
and 75400 s of exposure time in F606W and F814W, 
respectively. In the other two fields, we obtained 2 orbits 
through F606W, totaling 4800 s, and 4 orbits through 
F814W, totaling 9600 s. All images were calibrated in 
the HST pipeline with CALWP2 using OPUS version 
2006_6a for the 2007 data and 2008_5c for the 2009 data. 

To expand our radial coverage, we also reduced 2 fields 
in the o uter disk, previously s tudied and named S2 and 
S3 fsee iTikhonov et al.l 120031 for more details). These 
fields lie ~8' (7 kpc) SW of the nucleus and were taken 
as part of GO-5369. They contain 1200 sec of exposure 
in F606W each. Field S3 contains 4200 sec of exposure 
in F814W, while S2 contains only 2100 sec in F814W. 

2.2. Reduction 

The data reduction and photometry for the AN GST 
survey are fully described in lDalcanton et all (|2009|) . For 
convenience we provide a brief summary of the tech- 
niques here. 

The photometry was measured simultaneously for all 
of the objects in the unco mbined images u sing the soft- 
ware package HSTPHOT (jDolphinl |2000T ) . This pack- 
age is optimized for measuring photometry of stars on 
WFPC2 images using the well-characterized and stable 
point spread function (PSF) calculated with TinyTim. 6 
The software fits the PSF to all of the stars in each indi- 
vidual frame to find PSF magnitudes. It then determines 
and applies the aperture correction for each image using 
the most isolated stars, corrects for the charge transfer 
efficiency of the WFPC2 detector 7 , combines the results 
from the individual exposures, and converts the mea- 
sured count rates to the VEGAmag system. 

The HSTPHOT output was then filtered to only al- 
low objects classified as stars with signal-to-noise (to- 
tal counts from the star to total noise) >4 in both 
filters. The list was further culled using sharpness 
(\F606W sharp + F8UW sharp \ < 0.27) and crowding 
{F606W crowd + F8UW crowd < 0.7). The sharpness cut 
was chosen based on the distribution of values in the orig- 
inal catalog. The crowding parameter gives the difference 
between the magnitude of a star measured before and af- 
ter subtracting the neighboring stars in the image. When 
this value is large, it suggests that the star's photometry 
was significantly affected by crowding, and we therefore 
exclude it from our catalog. Quality cuts based on the x 
values were also considered, but they were rejected when 
a correlation was found between \ an( i the local back- 
ground. Our final star catalogs contained 40793, 22332, 
and 33365 stars for the deep, NE, and SW fields, respec- 
tively. The archival outer disk fields S2 and S3 are much 
more sparsely populated. We obtained reliable photome- 
try for only 468 and 692 stars for S2 and S3, respectively. 
The final CMDs for the fields are shown in Figures and 

13 

6 http:/ /www. stsci.edu/software/tinytim/ 

7 July, 2008 formulae 
http://purcell.as.arizona.edu/wfpc2_calib/2008_07_19.html 
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2.3. CMD Fitting 

We measured the star formation rate and metallic- 
ity as a functio n of stellar age using the software pack- 
age MATCH (jPolphinl 12003) . We fit the observed 
CMDs (with magnitude cuts set to limits provided in 
Table [B by popula ting the stellar evolution models of 
iGirardTet all rtgOOa with updates in lMarigo et~aT1l2008ft 
with a Salpeter (1955) initial mass function (IMF) for a 
grid of assumed distance and foreground extinction val- 
ues to allow for systematic differences in stellar evolu- 
tion models and/or systematic photometric errors. The 
choices of software and models used for the ANGST 
project are discusse d in detail in I Williams et all (|2009aj) 
and summarized in iDalcanton et al.l ( 2009TT 

The best fits provide the combination of ages and 
metallicities that are contained in the observed field. We 
attempted to fit the data with a spread in the model 
photometry along the reddening line to account for the 
effects of differential reddening. However, applying a 
spread in reddening of Ay — 0.5 to the models de- 
graded the quality of the CMD fits, showing that differen- 
tial reddening does not significantly affect our measure- 
ments in NGC 404. The data from the deep field were 
best fit by a single foreground reddening Ay=0.1±0.06 
and m - M =27.48±0.09 (see Figure HJ. This distance 
modulus is consistent with, but larger than, the value 
measured by the ANG ST survey (m - M =27.42±0.02 
IDalcanton et al.l l2009f > . The best fit values compensate 
for any systematic differences between the data and over- 
all model isochrones, whereas the survey value isolated 
the well-determined location of the tip of the red giant 
branch in order to measure the best distance. Since these 
distance and extinction values provided the best over- 
all fit of the models to the data, we performed our fits 
to the data of every region assuming these values. Our 
uncertainties in star formation rate account for changes 
in the SFH measured if the assumed value for the dis- 
tance modulus was ±0.15 mag away from the chosen 
value and if the extinction value was ±0.1 mag away 
from the chosen value. We note that in Figure [3] the ap- 
parent magnitude of the tip of the RGB appears fainter 
than in our more populous CMDs. Our fits to these 
CMDs (Figure [5]) were not significantly improved by al- 
lowing a greater distance modulus, indicating that the 
small numbers of stars in these fields cause the tip of the 
RGB to be under-populated. On the other hand, the 
SFHs of these fields could be significantly different than 
the inner fields. The shallow depth and low numbers of 
stars in these fields limits our ability to constrain the age. 
Our fits show only that 99% of the stars are older than 
1.6 Gyr. Thus, it is possible that most of these stars 
are only a few Gyr old, which could produce a fainter 
TRGB. Unfortunately, there is not enough data to con- 
strain whether the cause is undersampling or age. The 
CMD can be fitted equivalently well by either possibility. 

Systematic errors are determined by the MATCH pack- 
age by comparing the results of SFHs from fits to the 
data with different values for the distance and foreground 
reddening to the field. These errors are then added in 
quadrature to the random errors governed by our sam- 
pling of the CMD. The random errors are determined by 
randomly drawing from the observed CMD to produce 
Hess diagrams that vary due to the Poisson statistics of 



our photometric sample. By producing and fitting 100 of 
these Monte Carlo CMDs, we are able to determine the 
rms of the residuals between SFHs from these fits and 
those from the fits to the original data. The combined 
la error measurements therefore account for the uncer- 
tainties in the distance to the galaxy, the foreground red- 
dening, any systematic shifts between the model colors 
and magnitudes and our measured photometry, as well 
as the number of stars and features present in our CMD. 

Our Monte Carlo tests are also used to determine our 
time sensitivity (i.e., the de pendence of our SFH o n our 
choice of time binning; see William s et al.ll20 09b. for a 
full description of the technique). Briefly, we calculate 
the standard deviation of the maximum likelihood value 
from our 100 runs. We assume that any fit to the data 
more than one standard deviation away from the best 
fit is unacceptable. We then rerun our fits while sup- 
pressing star formation in various time bins. If the fit 
quality does not change significantly (by more than one 
standard deviation), we continue to expand the length of 
these removed time bins until the software can no longer 
find an acceptable fit. At this time resolution, we can be 
confident that our data provide meaningful constraints 
on the SFH. The final time bins are all sensitive enough 
that their removal from the SFH results in an unaccept- 
able CMD fit. 

2.4. Depth 

Photometric depth determines the precision with 
which we can recover the SFH of a region. The ef- 
fects of age and metallicity are more difficult to distin- 
guish with shallow photometry than with deep photom- 
etry. Our deepest photometry comes from our deep field 
and reaches the red clump in the least crowded region. 
This photometry therefore provides the most leverage for 
breaking the degeneracy between the age and metallicity 
of the old populations. In contrast, the data from our 
2 shallow fields provide the least of this leverage. How- 
ever, since the stellar populations should be well mixed 
at ages 3>1 Gyr, we used the metallicity distribution for 
the old stars (>2 Gyr) as determined from the fit to our 
full deep field to limit the range of allowed metallicities at 
each age in the fits to the shallower data. When the free 
parameters used to fit the shallower data were limited, 
we found the resulting age distribution of the ancient 
populations of the shallower fields to be consistent with 
those of the deep data, and the quality of the fit remained 
in the acceptable range (within 1 standard deviation of 
the value obtained when the full grid of free parameters 
was allowed). In what follows, the full- field SFHs for the 
shallower fields are the best fits possible with the restric- 
tion that the ancient population (>2 Gyr) contain only 
the metallicities at each age that contributed to the best 
fit of the deep full-field data. 

2.5. Field Division 

In addition to the full field CMDs, we divided our in- 
ner disk fields into 4 annuli to look for radial trends of 
the stellar populations and to help minimize differential 
crowding effects across each field. The annuli were chosen 
so that each region contained ^10000 stars in our deep 
field data. Their parameters of the annuli are provided 
in Table [TJ Figure [5] shows the CMDs of the different 
annuli in our deep field. These CMDs show clearly that 
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crowding in the inner regions is significantly worse than 
in the outer regions of our fields (as shown by the lim- 
ited depth of photometry in the inner fields) . In order to 
constrain the old populations of the 2 inner annuli, we 
limited the number of free parameters by allowing only 
the metallicities that provided the best fit to the outer- 
most annulus of the deep field for ages >2 Gyr. This 
restriction produced a acceptable fits to the CMDs while 
forcing the ancient population of the innermost regions 
to be of similar metallicity to that of the outermost re- 
gion. 

2.6. Spheroid Contamination 

Bag gett et al.l (|1998|) performed an early bulge-disk de- 
composition on NGC 404 from digitized photographic 
plate observations in the F-band, finding a bulge effec- 
tive radius of 63.8" and a disk scale length of 129.5". 
Their parameters indicate that the disk population is 
responsible for 50% of the galaxy light at —165" (near 
the inner boundary of our third annulus from the cen- 
ter) . Adopting this decomposition, our data are centered 
around the transition region between the domination of 
the bulge component and the domination of the disk com- 
ponent, suggesting our 2 outer annuli probe regions that 
are dominated by the disk population. However, more 
recently, the light profile of the bulge was measured with 
HST data, and the resulting effecti ve radius was found 
to be smaller f38". lSeth et~aLll2010f) . We performed our 
own fit to the light profile using the Spitzer 3.6/im im- 
age from the local volume legacy (LVL) survey 8 for the 
central portions and normalizing to the star counts from 
our fiSTdata. We found a similar effective radius (42"). 
Our fit requires that 75% of the stars belong to the disk 
at 165". These measurements suggest that the amount of 
sphero id contamination inferred from the Baggett et al. 
(1998) decomposition is an upper-limit, and all of our 
data are dominated by disk stars. 

3. RESULTS 

3.1. The Ancient Stellar Population of the Disk 

The SFHs for our 3 fields are shown in Figure (SFR 
versus time) and Figure [5] (the cumulative fractional 
star formation versus time). The SFH of the outer disk 
(archival fields S2 and S3) was not well-constrained due 
to the small number of stars. We note only that the weak 
constraints obtained from these fields (>98% of the stel- 
lar mass formed before 1.6 Gyr ago) were consistent with 
the age distributions of the stars in the outer disk being 
the same as those in our 3 inner disk fields. 

The stellar populations are clearly dominated by stars 
formed prior to 10 Gyr ago in all 3 fields. The cumu- 
lative age distribution of all of the fields is consistent 
with ^ 70% of the stellar mass formed by z ~ 2. Most 
of this constraint comes from the deep data, which pro- 
vided the metallicity constraints for the shallower fields 
and has the greatest fraction of stellar mass formed prior 
to 10 Gyr. This initial star formation was then followed 
by roughly continuous star formation at low levels to 
very recent times ( H 10 Myr ago). GALEX ultraviolet 
imaging is consistent with this low level of star forma- 
tion continuing through the past few Myr and possibly 

8 http://irsa.ipac.caltech.edu/data/SPITZER/LVL/ 



to the present day. The low recent star formation rate 
is consiste nt with the low su rface density of H I (~1.2 
Mq p C - 2 . ldel Rio et al.ll200l . 

Our results show that the stars in our disk fields are 
older than any other disks previously-studied in this 
manner; however, it is not clear if all parts of our fields 
are uniformly old. To address this question, we measured 
the SFH of our radial annuli. The cumulative SFHs of 
the stars in each of the 4 annuli are shown in Figures [9] 
and[TU] These SFHs do not generally constrain the older 
populations as well as the full-field fits because the CMDs 
contain fewer stars. However, they are all consistent with 
£ 70% of the stellar mass being formed by z ~ 2. The 
deepest photometry in our sample, which comes from the 
outermost annulus of the deep field, constrains Si 90% 
of the stellar mass to be older than 8 Gyr (z — 1). The 
result from the outermost annulus of the shallower NE 
field is consistent with that result but does not provide 
nearly as tight of a constraint (>70% of the stellar mass 
older than 2 Gyr) due to the lack of any of the red clump 
feature in the CMD. In the end, the SFHs of our differ- 
ent annuli allow the possibility that the inner annuli are a 
few Gyr younger on average than the outer annuli; how- 
ever, it is not clear whether this difference is due to bulge 
contamination, shallower photometry, or true differences 
in age. 

Because the full SFH measurements in our radial an- 
nuli did not show any conclusive differences between the 
stellar populations of the NGC 404 disk with radius, we 
attempted a more traditional method o f comparing the 
AGB/ RGB ratios in different regions. iTikhonov et al.l 
(2003) found that inside of r -150", the ratio of AGB 
to RGB stars decreases with radius, which interpret as 
evidence for most of the AGB population being associ- 
ated with the bulge. Our constraints on the detailed 
age distribution in the inner annuli allow the possibility 
that the inner regions are younger than the outer regions. 
Furthermore, our innermost annulus has a slightly bluer 
RGB and a higher AGB/RGB r atio than the other an - 
nuli studied, consistent with the ITikhonov et al.l (2003) 
results. Gaussian fits to a slice through the CMDs in 
Figure [6] at mpsiAW = 24.75 give a central RGB color of 
mF606W — = 0.95±0.01 in our innermost annu- 

lus, while those of the outer annuli and outer disk fields 
are all >0.98. The ratio of AGB stars (mpsuw < 23.5) 
to RGB stars (25.0 > m F 814W > 23.5) is 0.056±0.005 in 
the innermost annulus, while those of the outer annuli 
are all <0.05. This ratio was not reliably constrained in 
the outer disk fields (S2 and S3), due to the low numbers 
of stars. Overall, the results of our full CMD fitting, as 
well as direct measurements of o ur CMDs' features, are 
consistent with the suggestion in ITikhonov et ail (2003) 
that NGC 404 contains a bulge population that is on av- 
erage younger than the disk. However, with the current 
data we cannot distinguish between this possibility and 
the possibility of radial variations in the mean age of the 
disk itself. 

The metallicity history from our deep field is shown 
in Figure 1111 As with many other nearby galaxies, we 
find that NGC 404 enriched to [M/H] > -1 quickly and 
then appears to have undergone more modest subsequent 
chemical enrichment. The errors are very large on the 
metallicity of the young populations both because the 
number of young stars is low and because the prop- 
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erties of the upper main sequence is not strongly af- 
fected by metallicity. Therefore, while our metallicity 
history does not show any obvious discontinuities, if re- 
cently accreted gas was of different metallicity than the 
previously-existing gas disk, it is not clear we would have 
detected the difference with statistical significance. In- 
deed, the allowed metallicity range in the 30-300 Myr 
timcbin allows any metallicity [M/H]>-1.7. 

From these results we can assign the typical age and 
metallicity of the stars in these portions of NGC 404 as 
12±2 Gyr and -0.75±0.37. These values are consistent 
with those determined from fits to optical spectra of the 
bulge (>50% of the stars >5 Gyr old and [M/H] —0.4 
ISeth et ahl |2010| ) , indicating overall similarity between 
the bulge and disk populations in NGC 404, partially 
explaining the lack of radial trends in our SFHs. 

The typical age and metallicity of the stars in NGC 404 
can be comp ared to those o f other nearby galaxies, as 
was done in IWilliams et al.l ()200 9a). This comparison 
is shown in Figure [12] It shows the NGC 404 disk fields 
have the oldest typical age of any disk fields in the nearby 
sample; however, the age was determined closer to the 
galaxy center than for any of the other galaxies. Note 
that the 2 oldest populations, those of the Galactic thick 
disk and NGC 404, are those measured closest to the 
galaxy centers. 

3.2. The Young Stellar Population 

Although the star formation in NGC 404 over the past 
Gyr amounts to an insignificant fraction of the stellar 
mass of the disk, it is of interest for understanding any 
radial trends in star formation activity as well as the ori- 
gin of gas in the NGC 404 disk and its likely fate. We 
show the most recent Gyr of the SFHs of the full fields 
in Figure Q2] along with the weighted average of all star 
formation in all fields. We show these results at a higher 
time resolution (^200 Myr), since our time resolution 
is good at these ages due to the rapidly-evolving bright 
main- and blue He burning sequences, which contain in- 
formation for ages ^, 500 Myr. For the outer disk fields, 
we were not able to reliably determine ages from 400- 
1000 Myr, due to the small numbers of stars in these 
fields. We therefore limited the age range to just two 
time bins for those fields. 

There do not appear to be any radial trends nor any 
strong positional dependencies of the recent SFH in the 
our NGC 404 disk fields (see discussion below). How- 
ever, a distinct pattern emerges in the global SFHs of 
the full fields as well as all of the annuli beyond the in- 
nermost (and shallowest) one. The pattern is plotted in 
Figures [T3] and [14] which show all of the recent SFHs in 
color and their mean in black. The best fits to all fields 
suggest that star formation increased substantially ~500 
Myr ago, rising from consistent with zero to ~2xl0 -4 
M Q yr _1 kpc -2 by ^400 Myr ago, but then falling to 
a lower level by the present. Indeed, no star formation 
is required during the interval 1.3-0.3 Gyr ago to pro- 
duce statistically-acceptable fits to the CMD of the deep 
field; likewise no star formation is required during the 
interval 2-0.3 Gyr ago to fit the NE field and during the 
interval 1.3-0.2 Gyr ago to fit the SW field. Likewise, no 
star formation more recent than ^100 Myr is required 
to fit the CMDs of any of the fields either, showing that 
any very recent (<100 Myr) star formation has been at 



a very low level, consistent with the lack of any diffuse 
Ha emission (Figure [T5|). Thus all of our measurements 
(outside of the innermost annulus) require star formation 
in the past few hundred Myr. While our data allow the 
possibility that this star formation was confined to the 
past 300 Myr, the best fits to our data have this star for- 
mation confined mainly to a small event that occurred 
~600-200 Myr ago. 

The recent SFHs from the outer disk data (archival 
fields S2 and S3, Figure [To| are consistent with those ob- 
tained from the inner disk after scaling by stellar density. 
However, we were unable to reliably constrain the age 
distribution earlier than 400 Myr ago due to the shallow 
depth and small number of stars in thes e fields, which 
are lo cated ^2.3 scale lengths (r s =129.5 // iBaggett et al.l 
1998) farther out in the disk. The best fits to outer disk 
fields suggest an increase in star formation beginning at 
a similar time (~400 Myr ago) as seen in the inner disk 
fields, with an intensity that is a factor of 10 less than 
in the inner disk fields, consistent with a simple density 
scaling. 

To look for variations in the most recent ( ^ 160 Myr) 
star formation with position in our inner disk fields, we 
computed the mass of stars formed in all 11 of our defined 
regions (4 annuli each in the deep and NE fields, and 3 
annuli in the SW field) during the past 160 Myr according 
to our measured SFHs. This age range probes all O- and 
early B-type stars, and thus provides a good baseline 
for comparison with the FUV fluxes measured from the 
GALEX data. 

We plot the median radius of the stars in each of our 
regions vs. the mass of stars formed in the past 158 Myr 
in Figure [17] No radial trend is seen, suggesting that 
overall the recent SFH of the disk has not been a strong 
function of radius in the inner (<5 kpc) disk, and is more 
consistent with star formation pe rcolating stochasticall y 
across the face of the inner disk (jMcQuinn et al.1 [20091. 
These results suggest that there is little, if any, change 
in the young disk stellar populations with galactocentric 
distance in NGC 404 from 1 to 9 arcminutes, which covers 
~3.7 scale lengths of the disk and ~9 effective radii. 

We have further tested our ability to detect the star 
formation responsible for the flux in the Galex FUV 
image, which appears as discrete knots of star forma- 
tion in the NGC 404 disk (see Figure □]). We first 
measured the FUV luminosity of these knots in each 
of our defined regions ( see Figure [T]) using SExtractor 
(jBertin fc Arnoutsl Il996| ). We then subtracted the ex- 
pected foreground/background luminosity as determined 
from a control region 10 scale lengths from the galaxy 
center. We find that the FUV luminosity does not corre- 
late with galactocentric distance (see Figure [151) , which 
is consistent with our result that the amount of re- 
cent ( ^ 160 Myr) star formation does not correlate with 
galactocentric distance (other than simply scaling with 
density). On the other hand, there is a hint of a cor- 
relation between the FUV luminosity of the knots and 
the mass of stars formed in the past 160 Myr, shown in 
Figure [TH] A Spearman Rank test reveals that a corre- 
lation is present at only the 2a level (5% probability of 
no correlation). 

The fact that the correlation is weak may seem surpris- 
ing, but the reason for the weakness is apparent when 
one looks at the mean age of the stars younger than 
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160 Myr in each region. The colors of the error bars in 
Figure [19] indicate this mean age, with redder colors rep- 
resenting older ages. With this additional information, it 
is clear that the most FUV-luminous regions contain the 
youngest stars, but not the highest stellar masses. At 
a given age (similar color points), the FUV luminosity 
correlates better with the stellar mass of young stars. 

Such a scatter in age should be refected in the FUV- 
NUV colors of the regions. We measured these col- 
ors and found ±0.8 mag rms scatter. About 0.3 mag 
of scatter is expected from photometric errors as de- 
termined by photon statistics. Some of the additional 
scatter is due to the presence of foreground stars in 
the NUV. Varying extinction is not severe enough to 
explain any additional scatter. In the GALEX band- 
passes, A FUV / E( B _ V ) =8.376 and Anuv / E(b-v) =8.741 
(|Wvder et all l20§5h . or E {FUV _ NUV) /E B _ V =-0.365. 
Thus, the scatter in E(b-v) needed to explain the mea- 
sured scatter in FUV-NUV is ±2.2 mag. Such a large 
amount of differential extinction would be far more than 
we observe in our CMDs. While we cannot easily quan- 
tify the effects of foreground contamination, the scatter 
in FUV-NUV color supports the idea that a range of 
ages is in large part responsible for the weakness of the 
correlation. 

Furthermore, the distribution of upper-main sequence 
stars on the GALEX FUV image, shown in Figure [15l 
suggests that these young stars are more widely dis- 
tributed than the FUV emission. While all FUV- 
bright regions contain upper-main sequence stars, Not 
all upper-main sequence stars are confined within FUV- 
bright regions. The different distributions could be ex- 
plained by stars migrating away from thei r birth regions 
on ti mescales <160 Myr (consistent with lLada fc Ladal 
2003), and FUV emission being domi nated by stars 
youn ger than 100 Myr (consistent with fGogartc n et al.l 
2009). These results show that there is a complex rela- 
tionship between UV properties, age, and stellar mass in 
regions containing recent star formation. 

3.3. The Origin of the Gas Disk 

It is interesting that our SFHs show an increase in 
st ar formation ~0-5 M yr ago, in light of the suggestion 
of Idel Rio et al.l (|2004l) that the gas disk may have origi- 
nated in a merger event with a dlrr galaxy 0.5-1 Gyr ago. 
However, our full-field CMDs all require some star for- 
mation between 2 and 10 Gyr, as well as star formation 
between 0.5 and 0.1 Gyr ago, to produce an acceptable 
fit to the observed CMD. Failure to include stars in the 
2-10 Gyr age range produces a model CMD containing 
more blue AGB stars above the RGB tip than are seen 
in the data. These stars produce a quantitatively signif- 
icant decrease in the quality of the model fit. The lower 
limits on the star formation rates in the 2-10 Gyr period 
are comparable to the peak rates over the past 500 Myr. 

To investigate the origin of the gas in the disk, we 
studied the correlation between gas density and star for- 
mation in the disk of NGC 404 over its history. We 
assume that all of the stars that formed during a given 
time interval in our SFH must have been in the form of 
gas at the beginning of that time interval. When this 
gas is added to the gas seen at the present day, we can 
infer the gas density in the past, and then compare it 
to the subsequent star formation rate in our SFH. We 



can do this at increasing lookback times, giving us a way 
of probing the correlation between gas density and star 
formation rate at a range of epochs. 

We begin with the current H I c olumn density of th e 
disk from 2-5' (~1.5xl0 20 cm" 2 ; Idel Rio et alll200l . 
which corresponds to 1.2 M Q pc~ 2 . We then assume 
that Sh = Shi a nd adopt the conversion E gas = 1.45Sh 
(|Kennicuttl I1989D to account for metals and molecular 
gas. Finally, we assume that all stars formed from a 
previously-existing gas disk that was in place at the be- 
ginning of each time interval. Our method therefore re- 
quires that the gas disk had a higher surface density in 
the past in order to provide the material for all of the 
star formation that subsequently occurred. We further 
assume that stellar evolution quickly recycles 20% of the 
mass of stellar mass back to the gas phase. This fraction 
corresponds to all stars >8 M (assuming a Kroupa IMF 
integrated from 0.1 M to 100 M ), which return their 
mass to the ISM on timescales shorter than 100 Myr. 
The exact value of this fraction had minimal impact on 
our estimates. 

With the above assumptions, we were able to infer the 
gas density of the disk for several epochs in the history of 
NGC 404 and compare that density to our measured star 
formation rates. We were then able to test the plausibil- 
ity of the assumptions by comparing the results to the 
known co rrelation between gas density and star forma- 
tion rate (jBigiel et al.ll2008D . The results of this recon- 
struction of the history of the gas density and star for- 
mation rate in the NGC 404 inner disk are shown in Fig- 
ure [201 where diamonds show the gas density and mean 
star formation rates of the disk for each of 5 epochs (0.3, 
0.6, 0.9, 5, 12 Gyr ago) calculated from our measured 
SFHs. The Ssfr vs. S gas in the early epochs agrees with 
the overall trends measured by the H I Nearby Galaxy 
Survey (THINGS) (jBigiel et al.ll2008l ) under the assump- 
tion that all of the stars formed from a gas disk that was 
formed ~14 Gyr ago. However, at 0.9 Gyr, the star for- 
mation rate is inconsistent with the gas density under 
this "closed box" assumption. 

To test the possibility that a recent merger with a small 
gas-rich galaxy may have affected the evolution of the 
NGC 404 disk, we performed a reconstruction of the past 
gas density and star formation rate described above, but 
we altered the closed-box assumption. Instead, we as- 
sume that 10 8 M of gas was added to the system 0.6-0.9 
Gyr ago. The effect of this assumption is to make the gas 
density of the disk significantly lower ~1 Gyr ago than 
it is today, unlike the closed box assumption where the 
gas density strictly increases with lookback time. The 
triangle points with dotted errors in Figure [20J show the 
results of this reconstruction. These points, including the 
one at 0.9 Gyr, are all consistent with the known corre- 
lation between gas density and star formation rate. This 
consistency allows the possib ility that gas was r ecently 
accreted, as was suggested by Idel Ri'o et al.l (|2004D based 
on the warped nature of the outer gas disk. 

Together, these results favor a scenario where a gas 
disk was in place for many Gyr. This disk was of 
low surface density (<5 M Q pc~ 2 ) for most of the 
history of the disk and formed stars at a low rate 
(<10~ 3 M yr _1 kpc~ 2 ). This gas may have been com- 
pletely depleted by ~1 Gyr ago, when new gas entered 
the system, either from surrounding gas filaments, or 
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from a small, gas-rich, dwarf galaxy that merged to pro- 
duce the current gas disk. This event triggered a small 
amount of new star formation over the past ^500 Myr. 

It is still possible that the present gas disk is a relic of 
the gas that originally formed the disk and that there was 
no recent merger. If so, the low-density gas disk must 
have formed stars in short, weak, episodes that lasted 
~200-400 Myr. Such episodes, like the one we have re- 
solved in the most recent Gyr, form only <, 10 7 M Q of 
stellar mass in the inner disk (r < 5 kpc, assuming ax- 
isymmetry), which could reconcile the amount of stars 
formed in the past Gyr with the amount of gas avail- 
able in our closed-box assumption. Such star formation 
would be insignificant compared to the £ 2xl0 9 M of 
stellar mass formed in the inner disk prio r to 10 Gyr ag o 
and the 1.5 xlO 8 M Q of H I present (|del Rio et al.ll2004f ). 
Such a scenario is favored by the isolated environment 
of NGC 404 but is not consistent with the star forma- 
tion rates we meas ure from 0.6-1.0 Gyr. Furthermore, 
I del Rio et al.1 {2004) suggest the warped nature of the gas 
disk as additional evidence for the recent merger. While 
the warp is consistent with the merger scenario, warping 
can also occur from misalignment of an gular momentum 
during late gas accretion from filaments {Jiang fc Binnevl 
1999t lShen fc SellwoodlfeOOeft and therefore does not re- 
quire a catastrophic merger. 

4. CONCLUSIONS 

We have measured resolved stellar photometry for ~15 
arcmin 2 of the inner disk and ~10 arcmin 2 of the outer 
disk of the SO galaxy NGC 404. Detailed fitting of the re- 
sulting CMDs shows that the disk is dominated by stars 
older than ~10 Gyr, with all of our inner disk data con- 
straining the age of ~90% of the stellar mass to >8 Gyr 
and our best data constraining the age of 75% of the 
stellar mass to >10 Gyr. These results show that the 
difference between this disk and most later-type disks 
that have been studied in detail is seen in the old stellar 



population as well as the young. 

We found no trends between the age distributions of 
the young or old populations and distance from the 
galaxy center within th e obse rved regions. However, as 
scon bv lTikhonov et al.l {2003ft the mean RGB color and 
AGB /RGB ratio of the innermost regions are consistent 
with a bulge population in NGC 404 that is younger than 
the disk. 

Despite its low gas density, NGC 404 appears to have 
been forming stars throughout its history on timescales 
of Gyr, even in the inner disk (~1— 1.5'). The SFH and 
the current density of the gas disk are consistent with 
NGC 404 having evolved roughly in a closed box over 
much of its history, starting as a large disk of high sur- 
face density and approximately following the Schmidt 
law as its overall star formation rate and gas density 
decreased. Such long-lived disks have been observed in 
isolated galax ies before, as in the compact dwarf ADBS 
113845+2008 (jCannon et alJ2009f ). whose gas disk is also 
of low density. Perhaps these low density disks have very 
long lifetimes due to very low star formation efficiencies 
at these densities. Our SFHs indicate that star forma- 
tion in the low density relic disk may be episodic, going 
through minor episodes of star formation that last a few 
hundred Myr every Gyr or so. 

On the other hand, in NGC 404 there is some evidence 
that this passive evolution was slightly disrupted by an 
event 0.6-0.9 Gyr ago, possibly a small merger. Such a 
merger event would alleviate the necessity for very low 
star formation efficiency ~1 Gyr ago by allowing that the 
currently-observed gas disk was not fully in place at that 
time. 

Support for this work was provided by NASA through 
grants GO-10915 and GO-11719 from the Space Tele- 
scope Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Incor- 
porated, under NASA contract NAS5-26555. 
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TABLE 1 

Properties of the Designated Regions 



Region 


Rin (") 


Rout (") 


R ("\& 


F606W50 b 


F814VF 50 C 


deep annulus 1 


64.4 


124.0 


104.7 


26.6 


26.0 


deep annulus 2 


124.0 


160.0 


142.1 


27.3 


26.3 


deep annulus 3 


160.0 


190.0 


174.1 


27.8 


26.9 


deep annulus 4 


190.0 


268.4 


215.1 


28.0 


27.2 


NE annulus 1 


93.0 


124.0 


111.9 


26.8 


25.8 


NE annulus 2 


124.0 


160.0 


140.9 


26.9 


26.1 


NE annulus 3 


160.0 


190.0 


173.0 


27.0 


26.1 


NE annulus 4 


190.0 


250.5 


208.3 


27.2 


26.3 


SW annulus 1 


64.0 


124.0 


95.2 


26.6 


25.8 


SW annulus 2 


124.0 


160.0 


136.8 


26.9 


26.0 


SW annulus 3 


160.0 


190.0 


168.7 


26.9 


26.2 


S2 


404.8 


565.2 


468.3 


26.8 


26.1 


S3 


406.7 


560.6 


480.8 


26.9 


26.4 



a The median galactocentric distance of the stars in the the region. 
b The 50% completeness limit of the F606W data. 
c The 50% completeness limit of the F814W data. 
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Fig. 1. — The locations of our NGC 404 fields (§ 12.111 are shown on a 3-color image using GALEX far-UV (blue), GALEX near-UV 
(green), and H I (red) images. North is up. East is left. Fields are labeled with the names used in the text. The radial annuli used for the 
analysis (§ 12.5ft are also shown and labeled. 
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DEEP NE SW 




F606W-F814W(vegamag) F606W-F814W (vegamag) F606W-F814W (vegamag) 

Fig. 2. — The color-magnitude diagrams of our 3 WFPC2 fields (§ I2.2D . Contours denote the de nsity of points in are as of the plot that 
would otherwise be saturated. Overlaid on the deep field CMD are example isochrones from the Girardi et al. (2002) models shifted to 
■m - M = 27.48, Ay = 0.08 (from blue to red: [M/H]=-0.4 and log(age) = 7.3,7.6,8.0,8.3,8.6, followed by l og(age)=10.0 and [M/ H] = - 
1.3,-0.7,-0.4,-0.2,0.0, respectively). A gray dashed line marks the tip of the red giant branch as measured by (Dalcanton ct al. 2009). 




Fig. 3. — The color-magnitude diagrams of the 2 archival WFPC2 fields in the o uter disk (ij | 2,2H previou sly studied bv lTikhonov et al.l 
(2003). A gray dashed line marks the tip of the red giant branch as me asured by { Dalcanton ct al. 2003). We note th e TRGB appears 
slightly fainter in these fields, explaining why the distance measurement of Tikhonov ct al. (2003) was greater than that of Dalcanton ct al. 
(2009)'. 
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Fig. 4. — Our best-fit model CMD to the outermost annulus of the deep field data set f£j I2.3II , This fit is our most reliable as it comes 
from our deepest photometry. Upper Left: Our observed CMD cut at the magnitude limits used for the fit. Upper Right: The best-fit model 
CMD. Lower Left: The difference between the data and the model. Darker colors denote excess stars in the data. Whiter colors denote 
excess stars in the model. The range is from -19 (lightest) to +23 (darkest), corresponding to fractional errors of ±0.5. Lower Right: The 
statistical significance of the residuals shown in Lower Left. The range is -5cr (lightest) to +10<r (darkest). 
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Fig. 5. — Our best-fit model CMD to the deeper of the two outer fields (S3). Panels are the same as Figure [4] The range is from -5 
(lightest) to +8 (darkest), corresponding to fractional errors of ~^q q. Lower Right: The statistical significance of the residuals shown in 
Lower Left. The range is -1.5<r (lightest) to +4.2(7 (darkest). 
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Fig. 6. — The CMDs of the 4 annuli we defined for studying radial variations in the stellar populations (§[23}. The inset of the upper-left 
CMD shows the locations of the stars in each annulus on our deep WFPC2 field. Annuli were chosen to each contain ~10000 stars in the 
deep field. Radii of the boundaries are 64", 124.0", 160.0", 190.0", and 300.0". 
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Fig. 7. — The SFH of our 3 fields, normalized by total stellar mass in the field. This normalization produces units of fraction of the total 
field stellar mass produced each Gyr (Gyr -1 ). Blue, green, and red correspond to the results from the deep, SW, and NE fields respectively. 
The units can be converted to Mq yr -1 kpc~ 2 by multiplying the Gyr -1 value s by 0.038, 0.080, and 0.038, respectively. The SFHs of all 
fields agree and show most of the star formation occurring before 10 Gyr ago (§ 13.11 ). 
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Fig. 8. — The cumulative SFH of our 3 fields. The results from each individual field are shaded in light gray. Areas of overlap between 
the SFH of different fields are darker. The darkest gray denotes the SFH that agrees with all of the fields. The results from all fields overlap 
significantly and suggest a disk dominated by very old ( ^ fO Gyr) stars (£j 13. ft . 
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Fig. 9. — The cumulative star formation history of our 4 annuli in just the deep field. The results from each individual annulus are 
shaded in light gray. Areas of overlap between the SFH of different annuli are darker. The darkest gray denotes the SFH that agrees 
with all of the annuli. The results from all annuli overlap significantly, although the inner annuli do not provide as much constraint on 
the age distribution of stars older than 2 Gyr as the deeper outer annuli. The fact that the SFHs of all annu li ov erlap with that of the 
well-constrained outermost annulus suggests little radial variation of the stellar populations in the deep field (§ 13.11 1. 
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Fig. 10. — The cumulative star formation history of our 4 annuli. The results from each field are shaded in light gray. Areas of overlap 
between the SFH of the same annulus in different fields are darker. The darkest gray denotes the SFH that agrees with all of the data. Th e 
results from all fields overlap significantly, but the deep field constrains the stellar populations to be dominated by the oldest stars (§ 13.11 . 
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Fig. 11. — The chemical enrichment history of NGC 404 from our CMD fit to the deep field. Thick error bars mark the spread in 
metallicity, while thin error bars show the estimated uncertainty in the mean metallicity and spr ead. The ancient stars are dominated by 
metallicities >-1.0, and there is evidence of an increase at times more recent than ~1 Gyr (§ I3.lt . 
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Fi g. 12. — Metallic ity ran ge of the dominant s tellar populati ons in deep resolve d photometry for M31 d Brown et al.l I2OO60 . 
M33 lIBarker et al.ll2007h . M81 l|Williams et al.N2009ah . NGC 5128 l|Reikuba et al.l l2005h . NGC 3377 dHarris et aij 120071). NGC 2976 
HWilliams et al.H2009bl) . and NGC 404 (this work) are plotted along with that of the Milky Way thick disk llAllende Prieto et al.l 120061) 
again st several other properties (§ET}. Top: The metallicities as a function of the absolute K-band magnitude of the galaxy ( Skrutsk ie et alJ 
20061'). Th e solid and dashed lines show the luminosity-mctallicity relations determined by Trcmonti ct al. ( 2004, gas phase metallicity) and 



Mouhcinc (2006, stellar red peak metallicities of galaxy "halos"), resp ectively. These rel ations were converted from B-band and V-band 
to K-band using the Tully-Fisher calibrations of Vcrhcijcn (2001) and Sakai ct al. (2000). The Milky Way luminosity was calculated by 
applying V ro t= 220 km s -1 to the Tully-Fisher calibration of Vcrhciicn (2001). Middle: The populations' metallicities as a function of the 
radii at which they were sampled (normalized to the half-light radius of the galaxy). Bottom: The populations' metallicities as a function 
of their ages. 
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Fig. 13. — Colored error bars: The recent star formation histories of the three fields. Colors, normalization, and units are the same as 
in [7] Black histogram: Combina tion of all of the SFHs of the individual regions. The combination shows a significant increase in the star 
formation rate ~500 Myr ago 13.211 . 
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Fig. 14. — Colored error bars: The recent star formation histories of all of the annuli outside of the innermost, normalized by total 
stellar mass in the field. This normalization produces units of fraction of the total field stellar mass produced each Gyr (Gyr -1 ). Black 
histogram: Combinati on o f all of the SFHs of the individual regions. The combination shows a significant increase in the star formation 
rate ~500 Myr ago (§ I3.2[ l, The combined values can be converted from Gyr -1 to Mq yr — 1 kpc — 2 by multiplying the Gyr -1 values by 
0.033. 
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Fig. 15. — The locations of t he up per-main sequence stars from our HST images (cyan crosses) are shown on a composite GALEX 
far-UV (blue), H« (red) image (§[32}. The FUV-bright regions all contain these young stars, but the young stars are not confined to the 
FUV-bright regions. Also note the very low amount of Ho emission present. 
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Fig. 16. — Colored error bars: The recent star formation histories of the two outer disk fields, normalized by total stellar mass in the field. 
This normalization produces units of fraction of the total field stellar mass produced each Gyr (Gyr -1 ). Black histogram: Combination 
of the SFHs of the two outer disk fields. The SFH is consistent with an increase in star formation beginning ~400 Myr a go, a nd the star 
formation rate over the past 400 Myr is consistent with a simple scaling of the rate in the inner disk fields with density (§ 13.211 . 
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Fig. 17. — Surface dens ity o f stars formed in the past 158 Myr as a function of galactocentric distance in the NGC 404 disk. No trend 
with radius is apparent (§ 13.21 1. 
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Fig. 18. — Galex FUV surface brightness as a function of galactocentric distance in the NGC 404 disk. No trend with radius is apparent 
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Fig. 19. — Total stellar mass formed in the past 158 Myr as a function of total Galex FUV flux measured in the regions shown in FigurclTI 
defined by the WFPC2 field edges and radial annuli. Redder colors denote higher mean ages for the stars formed in the past 158 Myr. 
O nly a slightly significant trend is seen, partially due to the differences in the mean age of the young stellar mass in the different regions 
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Fig. 20. — Diamonds with solid errors: Gas surface density at previous epochs assuming all gas was in place 14 Gyr ago. Three epochs 
are plotted. Redder colors denote older ages (0.3,0.6,0.9,5,12 Gyr). Triangles with dotted errors: Gas surface density at previous epochs 
assuming a gas accretion event occurred 0.6-0.9 Gyr ago. Colors are the same ages as the diamonds. Sh aded Area: The regions of gas 
density - star formation rate density space covered by the relation measured by THINGS (Bigiel ct al. 2008). Note that the star formation 
rate was far too low for the assumed gas density 0.9 Gyr ago if we assume no gas has been accreted (light blue diamond; £j 13. 31 . 



